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DERIVE & Windows VISTA 


There were some rumors that DERIVE would not be supported by Windows Vista. We could 
give a VISTA-screenshot with DERIVE earlier. 

Then I received the message that DERIVE is running fine but the Online-Help cannot be 
opened because VISTA does not support the old HTML-format. 

Only two days later our DUG-member, DERIVE & IT-specialist Gunter Schodl told me that 
there is a Microsoft patch to overcome this incompatibility. He immediately downloaded the 
patch an installed it on his laptop computer. You can find the URL for the patch below. 

Heinz Hinkelmann asked for the URL and he sent the VISTA screenshot with DERIVE and 
the Online-Help open. Thanks to Gunter and Heinz. 



http : / /www .microsof t . com/ downloads/ details . aspx?displaylang=de& Family I D=6eb 
cfad9-d3f5-4365-8070-334cdl75d4bb 

X64 for 64 Bit platforms 

X86 for 32 Bit plattforms 

Just launch the file in Vista. Much luck, Gunter. 

Congratulations 

I have the pleasure to inform you that our vice president 
Barbel Barzel made a great jump in her career. Barbel - a 
DUG-member since its first days - made one of her first 
steps in using CAS in math education in Krems (“Taylor 
Series Expansions”) and has ever been a enthusiastic 
propagator of CAS. (Do you remember the “Nonhigh-flyers”, 

Barbel?) She is now appointed professor at the Institute of 
Mathematics and Informatics and its Didactics at the 
Pedagogical Highschool in Freiburg. We are sure that Baerbel will continue her work with all 
her energy and enthusiasm. We wish Baerbel and her family the best for the future and are 
looking forward to meeting Professor Barzel at any occasion. Much Luck! Josef and the 
DUG-family. 
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Liebe DUG-Mitglieder, 

endlich kommt der DNL (fast) zeitgerecht zu Ende 
des 1. Quartals zu Ihnen. Wir haben diesmal wieder 
ein sehr umfangreiches User Forum. Ich finde es 
besonders bemerkenswert, dass ein Problem aus 
dem alten DNL# 13 iiber dessen revidierte Neufas- 
sung so groBes Interesse gefunden hat. Inzwischen 
sind weitere Kommentare zu „R. Schorns Problem 
aus DNL# 13“ eingelangt. 

Mit MuPAD konnen wir einen neuen CAS-Gast 
begriiBen. Thomas Himmelbauer arbeitet schon seit 
langerer Zeit mit diesem Programm in seinem 
M-Unterricht. Er schickte uns eine Unterrichtsein- 
heit zum Heronschen Verfahren zum Wurzelziehen. 
Ich konnte es naturlich nicht lassen, seine schone 
Animation mit DERIVE und ITERATES nachzu- 
vollziehen. Der rekursive Ansatz eignet sich auch 
ausgezeichnet sowohl fur den Home Screen als 
Einstieg und fur den Sequence Modus als weiteren 
Ausbau auf den TI-Rechnem. Leider lasst sich hier 
die graphische Darstellung nicht so leicht - oder 
nur iiber das Programmieren - herstellen. Die Ani- 
mation mit DERIVE ist mir gelungen, weil wir von 
Peter Schofield ein schones Werkzeug zur Erzeu- 
gung von „Slideshows“ erhalten haben. 

Das ist aus einem Wunsch einer Kollegin entstan- 
den, die geme die schrittweise Approximation eines 
Kreises und damit eine Annaherung an k iiber ein- 
geschriebene regelmaBige Polygone visualisieren 
wollte. Der slider bar last sich im Zusammenhang 
mit dem VECTOR-Befehl oft nicht einsetzen, aber 
Peter hat diese Schwierigkeit mit einem einfachen 
Trick iiberwunden. Sehen Sie selbst, was man dar- 
aus machen kann. 

Ich habe in diesem Zusammenhang Riemannsum- 
men animiert. Im nachsten DNL wird Wolfgang 
Propper zeigen, wie dieses Thema mit TI-NSpire 
bearbeitet werden kann. 

Ich mochte Sie noch besonders auf die Grafiken 
von Pierre Charland hinweisen. Mir fallt die Aus- 
wahl aus seinem iiberreichen Angebot immer sehr 
schwer. Eine groBe deutsche Zeitung (FAZ?) wid- 
mete vor einiger Zeit eine Doppelseite spannenden 
impliziten 3D-Flachen. Ich werde diese Flachen in 
den nachsten DNLs mit verschiedenen Program- 
men erzeugen. 

Mit groBer Freude habe ich wieder einen Artikel 
von Don Phillips aufgenommen, der unsere reiche 
Sammlung von wertvollen statistischen Werkzeu- 
gen erweitert. Zum Schluss bedanke ich mich geme 
bei Johann Wiesenbauer, der - verlasslich wie 
immer -neue Titbits vorstellt, die sich wieder mit 
den Sudokus beschaftigen. Allen Autoren herzli- 
chen Dank! 


Dear DUG Members, 

I am very glad that I can deliver DNL#65 - nearly 

- in time. As you will notice there is a very ex- 
tended User Fomm. Especially remarkable is the 
fact that there was and still is so much interest in 
“R. Schorn's problem from DNL# 13” after putting 
the revised DNL# 13 on the website. 

We welcome MuPAD as our special CAS-guest. 
Thomas Himmelbauer has used this CAS for a 
while in teaching and he sent one of his numerous 
teaching units: Heron's algorithm for finding a 
square root. I couldn't resist to reproduce his ani- 
mation with DERIVE. This was only possible be- 
cause Peter Schofield provided a Slide Show Tool . 
The recursive procedure can be demonstrated on 
the Tl-calculators, too: as an introduction on the 
Home Screen and then using the Sequence Mode. 
Graphic representation of the process by an anima- 
tion is not so easy to perform and requires some 
programming. 

I was asked by a colleague how to apply the slider 
bar to visualise the approximation of a circle and 
subsequently of n by a sequence of inscribed regu- 
lar polygons. I knew that slider bars have problems 
applying on VECTOR- generated lists. But Peter's 
trick overcomes this difficulty. See, what can be 
achieved now using his versatile tool. As further 
examples I added the circumscribed regular poly- 
gons together with a visualisation of the approxi- 
mations of 7i and from another field of mathemat- 
ics, the Riemann sums. Wolfgang Propper will 
show how to animate Riemann sums with TI- 
NSpire in the next DNL. 

I'd like to call your attention to Pierrre Charland's 
graphics. It was not easy for me to make a choice of 
his rich supply of great pictures. Speaking about 
graphs, a well-known German newspaper (FAZ?) 
dedicated a double page to exciting implicit defined 
3D surfaces. I intend to present these surfaces in 
the next DNLs generated by various programs. 

It is a pleasure for me to have another contribution 
from Don Phillips, which extends our collection of 
valuable DERIVE statistics tools. We close this 
DNL again with Johann Wiesenbauer 's Titbits. He 

- reliable as ever - treated once more the Sudokus. 
Thank you all for your wonderful cooperation! 



Download all DNL-Derive- and Tl-files from 

http : / / www . aus t romath . at /dug/ 
http : / /www . bk- teachware . com/ main . asp?session=37 505 9 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wurmla 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: June 2007 

Deadline 15 May 2007 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Are all Bodies falling equally fast, J. Lechner 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Mathematics and Design, H. Weller, GER 

Diophantine Polynomials, D. E. McDougall, Canada 

Challenger Matrix Problems, G P Speck, New Zealand 

Graphics World, Currency Change, P. Charland, CAN 

Precise Recurring Decimal Notation, P. Schofield, UK 

Problems solved using the Tl-Nspire, K. Stulens, BEL 

Step function and Riemann Sums on the TI-NSpire, W. Propper, GER 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ 

Centroid of a Triangle - Recursively, D. Sjostrand, SWE 
BooleanPlots.mth, P. Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Baumann, GER; Keunecke, GER, 
and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Johann Wiesenbauer, Vienna Tue 06.02.2007 15:25 

Hi folks, 

Now that Josef Boehm is about to reissue the DNL #13, 1 had a look at R. Schom's problem on 
page 3, namely to compute the maximal coefficient of the polynomial 

(4x 3 + 3x 2 + 2x + 1) 20 . 

Well, that was almost 13 years ago and it took Derive 387.2s then to find the answer 
8842311087597693745, which turns out to be the coefficient of x A 40. 

Just to take into account the advances of both Derive and computers since then, I would like to in- 
crease the exponent to say 100, and pose this as a new challenge. (As there are vacations at the univer- 
sities right now, I thought, you might feel like a challenge!) In other words, what is the maximal coef- 
ficient in the expansion of 

(4x 3 + 3x 2 +2x + 1) 20 

and in which monomial does it occur? 

I for my part also got my teeth into this nice problem and just in case you want to compare with my 
solution (in Derive 6.10), you will find it in the attachment. 

Cheers, 

Johann 


OStD. R. Schorn, Kaufbeuren, Germany 

In the ’’Journal of Recreational Mathematics” (Volume 25, Number 2 - 1993) I found the following problem: 
2052. A Coefficient Problem by Charles Ashbacher, Cedar Rapids, Iowa. 

On page 473 of the unsolved problems section of Index to Mathematical Problems: 1989-1984 by Stanley Ra- 
binowitz, there is the problem: 

Menemui 5.2.1: Find the largest coefficient in the expansion of (1 + 2x + 3x^ + 4x^)20. 


This note started an interesting discussion between Johann and Valeriu Anisiu from Roma- 
nia. I'll skip Johann's first version and give the last one at the end of the discussion. 

The first reaction came from Spain: 

Ignacio Larrosa Canestro, A Coruna (Espafia) 

Hi Johann 

I cannot open the attached file, but I spend 2.17 to develop the polynomial, 0.454 to aprox, and a few 
seconds more to visually inspect it. 

Then the maximal coefficient is 

3982436672034106230248482366846039631376192603563093332660220654629851559396041954 
63780443664821745 ~= 3.982436672* 10 A 98 

and occurs in the term correspondent to x A 200. 

Best regards, 

Ignacio Larrosa Canestro, A Coruna (Espana) 
ilarrosa@mundo-r.com 
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Johann Wiesenbauer, Vienna 

Hi Ignacio, 

Yes, your result is correct. It was meant to be a programming challenge though and getting the final 
result somehow doesn't matter here as much as the total time needed for this. Sorry, if my post had 
been misleading in that respect! 

Then Valeriu joined the party: 


Valeriu Anisiu. Romania Tue 06.02.2007 19:57 

Hello Johann, hello everybody. My solution is slower but simpler: 

niaxterml(u , n, x, c_, v_) := 

Prog 

#1: v_ := TERHS(EXPAND(u A n), x) 

c_ := SUBST(v_, x, 1) 
v_j. POSITION ( MAX (c_) , c_) 


However, the result of 

3 2 

#2: raaxterml(4-x +3-x +2-X+1, 100) 

namely 


#3: 39824366720341062302484823668460396313761926035630933326602206546298515593960419546- 


3780443664821745- 


200 

x 


is obtained in < 2 seconds. 


Best regards, 
Valeriu 


Another polynomial specialist appeared: 


Danny Ross Lunsford 


#1: PROG 


Tue 06.02.2007 18:23 

' b 


3 2 100 

a := EXPAND((4*x + 3-x + 2-x + 1) ), b := TERMS(a) , c := VECTOR 




301 - i 


C x 




\ 


i , 1, 300 


, DISPLAY(MAX(c)) 


y y 

This takes a little over 3 seconds on my ThinkPad A31p. 
Answer 


39824366720341062302484823668460396313761926035630933326602206546298515593960419546- 

3780443664821745 


101st term 
-drl 



D 
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Again Johann and then two replies from Valeriu: 


Johann Wiesenbauer, Vienna Tue 06.02.2007 18:26 

Hello Valeriu, 

Great! Yes, as for brevity and elegance, your program is hard to surpass. 

In particular, it makes excellent use of some built-in Derive features. Even though, I hate to say it, but 
it is considerably slower (by a factor of about 100 on my machine) than my program. 

You know, I have been in the programming business for decades now, and nowadays I'm no longer the 
merciless performance hunter I used to be. If a program is so elegant and short like yours, I would 
certainly accept a lower performance to some degree. But here, as for me, this sacrifice of performance 
is too much though. (Just a hint: Try to increase the exponent n in this example up to say n=1000 and 
compare!) 

Cheers, 

Johann 


Valeriu Anisiu , Tue 06.02.2007 22:16 

Hello Johann, 

I was aware from the beginning that your program is much faster. 

But your program does not use Derive as a CAS; it manipulates lists instead of polynomials (except 
for a derivative, but this can be also avoided). 

In an ideal CAS (an world) the manipulation of polynomials should be very efficient. My program is 
waiting for that moment -:). 

Best reards, 

Valeriu 

Valeriu Anisiu Wed 07.02.2007 00:58 

Hello Johann, 

Another thing, 

Your program seems to give wrong answers and becomes much slower than mine when the degree of 
p is » 4. 

Try maxterm(x A 100 + 4x A 3 + 3x A 2 + 2x + 1, 5) or more dramatically 
maxterm(x A 1000 + 4x A 3 + 3x A 2 + 2x + 1, 5), versus maxterml (...). 

Cheers, 

Valeriu 

Johann ended the discussion and his file summarizes the process of finding the solution: 

Johann Wiesenbauer, Vienna Tue 06.02.2007 23:27 

Hello Valeriu, 

Yes, you have a point there. If the degree of the base polynomial is very large, a direct application of 
my method fails. I don't think that this case is very typical though. 


Cheers, 

Johann 
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Johann Wiesenbauer, Vienna Wed 07.02.2007 11:36 

Hello Valeriu, 

Based on your remarks, I made a revision of my program (cf. the attachment), which looks much more 
streamlined now. Even though, your objections remains valid. In particular, as to maxterm() I have to 
add certain requirements regarding the base polynomial, otherwise it won't work correctly. (I'm sure, 
one could easily fix this problem, but I don't have the energy right now. One should remember that I 
wrote that routine to solve my original challenge efficiently, which it does, indeed!) 


At any rate, many thanks for your valuable remarks! 


Cheers, Johann 

polypowers . dfw 


Yet another way of computing polynomial powers 

(by Johann Wiesenbauer, Vienna University of Technology, revised on February 7) 

If u is any nonconstant polynomial in x, then the following routine computes the coefficient of 
x^k of the power y = u^n for any integer n > 1 by solving the differential equation u y' = n u' y 
for y in an iterative way. 

coefffu, n, k, a_, b_, c_, i_ := 1, m_) := 

Prog 

m_ := POLY_DEGREE(u) 

a_ := ADJOIN (SUBSTfu, x, 0) A n, VECTORfQ, j_, m_ - 1)) 

b_ := VECTOR (P0LY_C0EFF(u , x, k_) , k_, Q, tO 

c_ := VECTOR + j_ - + 1), j_, 1, m_) 

# 1 : Loop 

If i_ > k 

RETURN FIRST (a_) 

a_ := ADJOIN (a_-c_/( i_. FIRST OO), DELETE (a_, -1)) 

c REST (b_) 

i_ :+ 1 

3 2 

#2: coeff(4-x + 3-x + 2-x + 1, 20, 40) = SS423110S7597693745 

In particular, we can use it in a slightly modified form to compute the monomial in u^n with the 
maximal coefficient under the assumption that all coefficients of u are nonnegative and nonde- 
creasing w.r.t. their index (see below). This solves a problem posed by R. Schorn in DNL #13, p. 3 
in a very efficient way from an algorithmic point of view. (Note the computation time of 0.000s 

r\y\ m\/ 

maxtermfu, n, a_, b_, c_, i_ := 1, m_) := 

Prog 

m_ := P0LY_DEGREE(u) 

a_ := ADJOIN (SUBST (u , x, 0) A n, VECTOR (0, j_, m_ - 1)) 

b_ := VECTOR CP0LY_C0EFFCu, x, k_) , k_, 0, m_) 

c_ := VECTOR ( (j_ ■ n + j_ - l)-b_j,Cj_ + 1), j_, 1, m_) 

#3: Loop 

m_ := a_- c_/(i_- FIRST (b_)) 

If m_< FIRST CO 

RETURN FIRSTCO->! A Ci_ - 1) 
a_ := ADJOIN Cm_, DELETE (a_, -1)) 
c_ REST (b_) 

i_ :+ 1 

3 2 40 

max term (4 -x + 3-x + 2-x + 1, 20) = SS423110S7597693745 -x 


# 4 : 
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The actual superiority of the algorithm above only shows though, if n is much larger, say 100, as 
in my challenge on the Derive-forum. 

3 2 

#5: maxterm(4-x + 3-x + 2-x + 1, 100) 

#6: 398243667203410623024848236684603963 1376192603 56309333266022065462985 15 59396041954~ 

200 

63780443664821745 • X 

As for this ongoing competition on the Derive-forum, the best solution so far was contributed 
by Valeriu Anisiu. It is very short and elegant and looks like this: 

maxtermlCu, n, x, c_, v_") := 

Prog 

#7: v_ := TERMS (EXP AND(u A n), x) 

c_ := SUBST (v_, x, 1) 
v_| POSITION CHAXCO, c_) 

Alas, as for the performance of the problem above, it is not a serious competitor being slower 
by a factor ~ 100. 

3 2 

#8: maxterml(4*x + 3-x + 2-x + 1, 100) 

#9: 3982436672034106230248482366846039631376192603563093332660220654629851559396041954- 

200 

63780443664821745 -X 

On the other hand, in a subsequent post Valeriu wrote: 


Another thing, 

Your program seems to give wrong answers and becomes much slower than 
mine when the degree of p is » 4. 

Try maxterm(x A 100 + 4x A 3 + 3x A 2 + 2x + 1 , 5) 

or more dramatically maxterm(x A 1000 + 4x A 3 + 3x A 2 + 2x + 1, 5), versus maxterml (...). 

Yes, as I already wrote in my answer, he has a point there. (Compare the run times below!) 


#10: 

100 3 2 

maxtermfx + 4-x +3-x + 2«x + 1, 

10 

5) = 17203- x 

#11: 

1000 3 2 

maxtermfx + 4-x +3-x + 2«x + 1, 

10 

5) = 17203- x 

#12: 

100 3 2 

maxtermlfx +4-x +3-x + 2«x + 1, 

10 

5) = 17203- x 

#13: 

1000 3 2 

maxtermlfx +4-x +3-x +2-X + 1 

10 

, 5) = 17203- x 


For bigger exponents the results become even wrong, because the requirement that the coeffi- 
cients of the base polynomial are nondecreasing w.r.t. the index is not fulfilled here. 


p 
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1000 5 2 20 

#14: maxterm(x + 4-x + 5-x + 2-x + 1, 10) = 1259550265 -x 

1000 5 2 101S 

#15: maxterml(x + 4-x + 5-x + 2-x + 1, 10) = 1505775150-x 

1000 5 2 

#16: coefffx + 4-x + 5-x + 2-x + 1, 10, 1018) = 1505775150 

Hence, when applying my routines, one should be aware of those restrictions. Many thanks to 
Valeriu for pointing out them! 


The problem was announced as unsolved in the 1980's. I don't believe that it was the inten- 
tion to have a computer solution. I wonder if there is a way to find this number by reasoning 
applying or developing something like a “super binomial formula”? Josef 


The Klein Bottle and its Normals 

Lester Anderson 
Hello 

I have a set of parametric equations defining a 3D surface, in this example a "Klein bottle", as a non- 
orientable surface. DFW file attached. 

k : =4 (1-cos (u) 12 ) 

For u=0-pi 

x=6cos (u) (1+sin (u) ) +kcos (u) cos (v) 
y=16sin (u) (1+sin (u) ) +ksin (u) cos (u) 
z=ksin (v) 

For u=pi-2pi 

x=6cos (u) (1+sin (u) ) +kcos (u) cos (v) 
y=l 6sin (u) 
z=ksin (v) 

v= -pi to pi 

What I would like to do is create a demonstration plot that shows normals to the surface (say at the 
mesh nodes for z=0, panels=20), to illustrate the nature of the ’’one-sided” surface i.e. no volume in 
this case. I am not too sure how to proceed on this other than setting a plane to z=0 as a starting point. 

Ideally I am looking to do a cutaway view. In the above case setting v=-pi to 0 for illustration. 

This is quite a complicated modelling task as not only do I need to determine the normals to the sur- 
face at a given z value, but also to make the normal vectors a small enough length to be instructive. 
Hopefully someone can provide some guidance on this issue as I would be looking to apply the tech- 
nique to other surfaces. 

The DFW file is in version 5 format, but works fine in version 6. Any advice about how to proceed 
will be welcome. 


Thanks 

Lester 


D 


-N-L#65 


Derive- and CAS-TI-User Forum 


P 9 


Dear Lester, 

here is bottle#2 + normals. 

Best regards 
Josef 

The plot is embedded - You can reanimate the Bottle, Cheers! 


# 1 : 


" 

r u ” 


' u ' 

X 



r u " 

a + COS 


SIN(v) - SIN 


•SIN(2*v) 

■COS(u), 

a + COS 


. \ 

v 2 , 


V 2 , 

j 



v 2 > 

( u ’ 




u ^ 

/• 

u 



SIN(v) - 


SIN 


■SIN(2-v) 


■SIN(u) , SIN 


■SIN(v) + COS 


0 <= u <= 7PI/4; 0 <= v <= 2PI; a = 3 


working with slider bars for a and b: 


#2: r := 4. 


1 - 


COS(u) 


■SIN(2-v) 



#3: [a • COS(u) • (1 + SIN(u)) + r • COS(u) -COS(v) , b-SIN(u) + r ■ SIN(u) -COS(v) , r-SIN(v)] 

#4: [a.COS(u).(l + SIN(u)) + r-COS(v + w), b.SIN(u), r.SINCv)] 



Lester Anderson's definition: 

#5: [ 6 • COS ( u ) • ( 1 + SIN(u)) + r -COS(u) -COS(v) , 16-SIN(u) + r.SIN(u).COS(v), r-SIN(v)] 

#6: [ 6 • COS ( u ) • ( 1 + SIN(u)) + r -COS(u) -COS(v) , 16-SIN(u), r-SIN(v)] 


Parameter curve u = PI/4 




pio 
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#7: 


#8: 


r 2 

LC0S(u)-(6-SIN(u) + 4-COS(v) + 6) - 2-COS(u) .COS(v), SIN(u)-(4-COS(v) + 16) 

2-SIN(u)-COS(u)-COS(v), 4-SIN(v) - 2 -COS(u) -SIN(v)] 

[ 2 

LC0S(u)-(6-SIN(u) + 4-COS(v) + 6) - 2-COS(u) .COS(v), 16-SIN(u), 4-SIN(v) - 

,] 


2 -COS(u) 'SIN(v). 


#9: 

c° s f— j 

• [e-sinf— 

-1 + 4-COS(v) + 6| - 2 • COS [ 

2 r n > 

•COS(v), SIN 


l 4 J 

l l 4 

J J l 4 J 

l 4 J 


•COS(v) + 


16) - 2 -SIN 


f n ) f ji ) 

-COS 

4 J { 4 J 


COS(v) , 4-SIN(v) - 2 -COS 


(-) 

l 4 J 


•SIN(v) 


#10: [(2-V2 - l)-COS(v) + 3-J2 + 3, (2-J2 - l)-COS(v) + 8-^2, (4 - 72)-SIN(v)] 

#11: par_c :i [(2-^2 - l)-COS(v) + 3-J2 + 3, (2-J2 - l)-COS(v) + 8-J2, (4 - /2)-SIN(v)] 

to produce a thick line: 

( 71 


#12: VECTOR [par_c], v, —77 , 77 , 


100 J 



finding the normal vector with respect to the surface 


#13: 


fd 

[6 * COS(u) ■ (1 + SIN(u)) + r * COS(u) ■ COS(v) , 16.SIN(u) + r • SIN(u) - COS(v) , 
^du 


r-SIN(v)] 

j 


d 

X — [6*C0S(u) • (1 + SIN(u)) + r -COS(u) -COS(v) , 

dv 


16 • SIN(u) + 


#14: 


r-SIN(u)-COSCv), r-SIN(v)] 

3 2 2 2 2 2 

8-COS(u) -COS(v) - COS(u) -(16-COS(v) + 32-COS(v)) - COS(u) • (4-SIN(u) -SIN(v) 

2 2 3 

- 12 -COS(v) - 64-C0S(v)) + 8-SIN(u) , 24-COS(u) -COS(v) + 
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Set again u = PI/4 





3 2 

' JT ' 

#15: 

8* COS 


•COS(v) - COS 




V 4 , 


v 4 > 


2 2 

- (16 • COS(v) + 32 * COS(v)) - 



f 71 ' 



r n ^ 

2 2 2 ^ 


r ^ 

COS 


* 

4- SIN 


•SIN(v) - 12 -COS (v) - 64 • COS(v) 

+ 8 * SIN 



V 4 > 


V 

v 4 > 

J 


k 4 > 


, [, 


2 2 
#16: L(9-72 - 8)-C0S(v) + (32-72 - 16)-COS(v) -72+4, (9-72 - 8)-COS(v) + (12-72 


6) • COS(v) +72-4, (18 - 8-72) • SIN(v) -COS(v) + (44 - 11-72) -SIN(v)] 
nv ist the normal unit vector of the surface for the points of the parameter curve from above: 


#17: nv := 


2 2 ~ 
(9-72 - 8) -COS(v) + (32-72 - 16)-COS(v) -72+4, (9-72 - 8)-COS(v) + (12-7- 


(9-72 - 8)-C0S(v) + (32-72 - 16)-COS(v) -72+4, (9-72 - 8)-COS(v) + (12-7- 

2 

(4-72 - 9) -COS(v) + 8- (72 - 4)-COS(v) - 2-72 + 1 


7((451 - 244-72) -COS(v) + 4-(222 - 115-72)-COS(v) - 492-72 + 1107) 
2 

(4-72 - 9) -COS(v) + 3 -(72 - 4)-C0S(v) +2-72-1 


2 

7((451 - 244-72) -COS(v) + 4- (222 - 115 -72) -COS(v) - 492-72 + 1107) 


SIN (v) • (11- (2 -72 - 1) - (8 - 9-72) -COS(v)) 


2 

7((451 - 244-72) -COS(v) + 4- (222 - 115 -72) -COS(v) - 492-72 + 1107) 


This creates a family of normal lines setting t from 0 to 3 after issuing Insert > Plot in the 3D PLot window 

n 1 

par_c + t-nv, v, -n, ji, 

6 J 



#19: VECTOR 


v. 



For the second calculation presenting another form of the Kleinbottle together with two pa- 
rameter curves and their normals resulting in the graph above example follow the file 
kleinbottle. dfw. Josef 


Lester Anderson 

Dear Josef 

Thanks for the help on the problem, it looks really good. It's good to have a starting point like this for 
topological studies; shows just what can be achieved with Derive! 


Peter Schofield 

Hello Joseph, 

I have just managed to download and have a good look at DNL#64, and I must thank you for your 
really nice presentation and selection of graphics examples using my “ImplicitPts” and “ImplicitDots” 
Derive files (DNL#64 pp. 33-40). 

I note that you have found some more interesting examples- 1 particularly liked the “Window of 
Viviani” formed by the intersection of a sphere and cylinder. I also found Tania Kollar’s article in 
DNL#63 (pp. 6-14) very fascinating. 

In 3D, using Derive Boolean expressions, it is also possible to generate “clouds” of points and 
dots using the attached short file: BooleanPlots.mth. Applications are limited by the size of the list of 
points (dots), but the file illustrates some of the possibilities. I’m sure that, given the opportunity, De- 
rive would have developed further graphic facilities in this direction. 

Many Thanks! 

Peter 


BooleanPlots.mth will be presented in one of the next DNLs, Josef 
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Plot Range for Parametric Plots 

Walter Wegscheider, Austria 

Dear Josef, 

this is a DERIVE related question: is it possible to influence the range of the parameter for parametric 
plots by a command. 

Example: vertical line for x = 2: g ( t ) : = [ 2 , t ] 


If I want to have 0 < t < 8 it is no problem to indicate the range in the 2D-plot dialogue. Unfortunately 
there is no command like g ( t ) :=[2,CHI(0,t,8)*t] or similar. Is there an alterantive via the 
status variables or . . . 


Dear Walter, 

Basically there is no way to pass the dialogue box. I offer a little trick (#17) which works prop- 
erly. Please inspect the attached file. 


*1: (TABLE([2 , t], t, 0, S, 0.01)) U [2, 3] 


#2 


: ^TABLE^t 2 ■ COS (t) + l, i-SIN(t) +4], t, s s G.0ljju[2, 3] 


Applying Point Size Medium or Large you can get thick graphs. 



pai'aiaf (pf , p_, start, end, step := 0.01) := (TABLE(Tpf , pf 1 , p_, start, end, step))||[2, J] 

U: 1 1 2l 

t4: paraiaf ([2 -SINH(ii) , - COSH(ii)], 11 , 0, 2) 

*5: paranf ([2 -SINH(u) , - COSH(ii)], 11 , -2, 2, 0.005) 

If necessary one has to change the step size. In the “standard plot window” step = 0.01 is 
sufficiently small. Don't forget to activate Approximate before plotting in the Plot window. 
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Finally I come back to the initial plots. #6 gives the plot of the vertical segment together with 
the segment of the ellipse. 


# 6 : 


f n 

parainf ([2 , t], t, 0, S) , parauf I [2 -COSCs) + 3, 3-SIN(s) +4], s, 5 


3 'J7 


4 


i 


Walter Wegscheider 

Dear Josef, , 

many thanks - your “trick 17“ means: program your personal plot command ©. 

The internal plot procedure works similar (including an extra check of the plot range considering the 
pixels of the screen). 

We have discussed this question several times and have the strong opinion, that both alternatives 
should be provided, which are in particular: 

a) the dialogue principle with intelligent default settings (like in DERIVE) which enables a first 
“beautiful 44 result without special syntax knowledge. 

b) the command principle (eg MuPAD) - each command must be entered with all possible parame- 
ters to obtain a reasonable result. This requires knowledge of the syntax but has the advantage 
that all settings can be saved in a document and can also be used in programs (and macros) be- 
cause all details can be addressed directly. 

Thanks again and best regards, 

Walter 

Annotations! 

James Gordon 

I have Derive 6. 

My expression #4 is sqrt(24). When I highlight #4 and execute the Command Toolbar's Approximate 
button, I get in the Status Bar's second pane "Approx(user)" rather than "Approx(#4). My settings are 
factory default. How can I get the Status Bar's second pane to read "Approx(#n), when n represents 
any item number. 

Thank you. 

Albert D. Rich 

Hello James, 

In order to get "Approx(#4)" use the Approximate button usually near the top of the Derive 
algebra window instead of the Approximate button to the left of the expression entry line. 
Hope this helps. 

Aloha, 

Albert D. Rich, Co-author of Derive 

Calculus Made Easy 

CME v.9.0 is available for all TI-CAS-Calculators 

http://www.ti89.com 

Our friend Nils Hahnfeld released a new version with many new and exciting features. More 
about CME v.9.0 in the next DNL. 


D 


-N-L#65 


Derive- and CAS-TI-User Forum 


pl5 


Francisco Marcelo Fernandez 

Dear Derivers: 

I am calculating Hankel determinants with the simple function 

dete(n, d) := DET(VECTOR(VECTOR(APPEND("c", i + j + d + 1), j, 0, n - 1), i, 0, n - 1)) 

where n is the determinant dimension, and d shifts the elements subscript. I obtain all the determinants 
with n = 2,.. .,7 easily and quickly, with (for example) d= 0. But the case dete (8,0) takes too long 
and consumes too much memory. Does anybody know why it happens?. I find exactly the same prob- 
lem with Derive versions 5 and 6. 

Thanks, 

Marcelo 


#1: dete(n , d) := DET(VECTOR(VECTOR(APPEND(c , i + j + d + 1), j , 0, n - 1), i, 0, n - 1)) 

2 2 3 

#2: dete(3, 0) = cl*(c3*c5 - c4 ) - c2 -c5 + 2-c2-c3-c4 - c3 

2 2 3 2 2 

#3: dete(4, 1) = c2 . (c4* (c6-c8 - c7 ) - c5 -c8 + 2-c5.c6.c7 - c6 ) + c3 .(c7 - c6*c8) 

2 3 2 2 

+ c3 • (2 «c4. (c5 *c8 - c6*c7) - 2-c5.(c5.c7 - c6 )) - c4 -c8 + c4 .(2-c5-c7 + c6 ) - 

2 4 

3 • c4 • c5 *c6 + c5 


A Normal Distribution Standard Problem 


Josef Lechner 

Sehr geehrter Herr Dr. Lechner, 

mit groftem Interesse habe ich Ihren Artikel "Standardisierung der Normalverteilung - ein Anachronis- 
mus?!" in den TI-Nachrichtenl/03 gelesen und stimme Ihren Urteilen zu. Bei meinem Unterricht zur 
Normalverteilung habe ich mich durch ihren Artikel anregen lassen und zwarauch die Gauftsche 
Dichtefunktion und ihre Beziehung zu anderen Dichtefunktionen normalverteilter Zufallsgroften be- 
handelt, fur Berechnungen wollte ich aber die Dichtefunktion benutzen, so wie Sie sie auch in Ihrem 
Artikel angeben. Die Berechnungen waren zunachst problemlos moglich bis wir auf folgende ganz 
normal und harmlos aussehende Aufgabe stieften: 

Eine Maschine stellt Nagel her. Die Lange der Nagel ist normalverteilt mit dem Erwartungswert p = 
8,00 cm und der Standardabweichung o = 0,15 cm. Bei wie viel Prozent der Nagel weicht die Lange 
hochstens urn e = 0,20 cm vom Erwartungswert p ab? 


Josef Lechner forwarded the following mail which he received from Dr. Weifi, a German Voyage User. 
Dr. Weifl refers to Josef's paper about “Standardisation of the Normal Distribution - An Anachro- 
nism?!”. Everything was ok, but one very common problem was a “problem” for the Voyage 200. 


A machine produces nails. The legth of the nails follows a normal distribution with mean p = 8,00 cm 
and standard deviation o = 0,15 cm. What percentage of the nails lies between p ±0,20 cm? 


As you can see, the calculator doesn 't deliver any result 
and shows a system message 

Overflow replaced by co or -oo. 

There were nor problems when working with the same 
data in mm or taking other values! 






* phi (u ? (j ? x) 


Done 


■J®|gPhi(3, .15,x)(ix 



WflKnin3: DutKflOW by to OK -to 


undef 
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Sie sehen: DerTi-92 meldete: Overflow replaced by unendlich or-unendlich und zeigte kein Ergebnis 
an. 

Der Rechner lieferte problemlos das richtige Ergebnis, wenn man die Zentimeterangaben aus der 
Aufgabe in mm umwandelte und ihn dann rechnen lieft. Ich habe keine Erklarung fur dieses Phano- 
men gefunden. 1st Ihnen im Unterricht mal Ahnliches passiert und/oder konnen Sie eventuell den Ef- 
fekt erklaren? Uber eine Antwort wurde ich mich sehr freuen. 

DNL: 

This is very interesting. Obviously the Tl does not like to accept this special value for the standard 
deviation (don't ask me why). Dr WeifJ (and we, of course) can overcome this bug by applying nu- 
merical integration - using nlnt. 

See the following Tl-screens. 

Interestingly enough I am lucky to produce the nice 457 digit number and infinite, but without an error 
message - so the result should be correct - or not? - ©. (It's only “questionable accuracy”). 


F£T Y FI:t 1 

t Algebra Calc 

f F4T 1 
Other 

r ?-■ i 

PrgnlO 

f FfiT 1 
Clean Up 

'“"I 


+ phitiu, x) 


Done 


■ nlnt(phi<8, . 15, x), x, 7.8, 8.2) .3176 

■J ® -2 phiCS, . 15, x!)cix 5.7750e382 

■j!' 2 phiCS, . 15, x!)cix 5.7750E382 

f <phi<8 --16-x>-x-7-8-8> 

CRYFTD FiliD EKflCT FUNC 1F/30 


[^7*01 

t i — 

fy. J 

F£T 

Algebra 

7 

f F ^ T 1 

Calc 

f F4t 1 
Other 

r ff i 

PrgnlO 

f FfiT 1 
Clean Up 



I 7.3 

[?.e 

i? 


phi(8 ? . 15 ? x)dx 
phi(S ? . 13 ? x)dx 
phi(S ? . 13 ? x)dx 


5. 7753 e3S2 
3. 3732 e457 


X<phi<8, .16,x>,x,7,9> 




Solids of Revolutions and Infinite Series 

Nils Hahnfeld, Virgin Islands 

(1) Short question: how can I show my students on the TI-89 a 3D-plot which results in a rotation 
of sqrt(x) between [0, 4] around the x-axis? 

DNL: There were several questions about representing solids of revolution with DERIVE 
and with the Tl-devices in the past. (See DNL#35). It is no problem to produce very pretty 
and informative plots applying the parameter form and then it doesn't matter which rotation 
axis applies. This is difficult on the TIs because we can only plot functions of form z = f(x,y). 
So we have to distinguish between the top half and the bottom half of the surface and we 
can plot only one at a time. 


So let's try to do our best: 

We rotate y = f(x) around the x-axis. 

The sketch explains the procedure: This is an 
intersection perpendicular to the rotation axis 
(= x-axis) at position x. The intersection is a 
circle with radius f(x). Point (x, y, V((f(x) 2 -y 2 )) 
lies on the top half of the surface and point 
(x, y, -V((f(x) 2 ->^)) on the bottom half. 



Take the TI-89 (or the TI-92 or the Voyage 200 ... or DERIVE): 

First of all I produce the paraboloid by parameter form with DERIVE: 
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left: parameter representation 
#1: [t, ^t-ODSCs), ^t-SINCs)] 


right: Two-functions-representation 


#2: 


2 2 

_IFCx < 4, ,/(x - y IF(X < 4, - J Cx - y ))_ 



The surface on the TI-89 screen: 



f Fl-r T F£t T F3 T FH T FE*T Ffi*TF7-*#i:'| 
|T-Mjls|Z-M:-rM|TKO-:h:|F;43i-OFh|Moi:h|DKOw|Firi|:-: | 






MAIN RAP AUTD 


3P 



Two More examples (V 200): Rotation of the sine wave 



J'LDTS 


-■zl=J(sin(;x)) 2 - y 2 


z2= -J(sin(x)) 2 -y 2 

z3= 

z4=l 

z5= 

zb= 

z7= 

zS= 


z4<x,y>= 


MAIN RAD AUTD 3D 


pi=q 

f F£T Y 1 

Zoon| | 

eye9= 



eye^=2G. S9 
xnin= "1 . 
xnax=4. 
xgrid=4G. 
ynin= "1 . 2 
ynax=l . 2 
ygrid=4G. 
znin= "1 . 2 
znax=l . 2 
ncontour=5 

i 


MAIN RAD AUTD 3D 





MAIN RAD AUTD 3D 


z2 shows the other half of the surface. Finally the rotation of f(x) = 2 - x A 2/2 for -3 < x < 3 




zl = .|( £in(x:)) 2 -y 2 

z2= -Ji;sini:^) 2 -SJ 2 


^z4 ; 

z5= 


z4<X . V>= zSCx^V^ 


eyeGH 

eye+= 

ey&T=22. 91 

xnin= "3. 

xnax=3. 

xgrid=4G. 

ynin= - 3. 

ynax=3. 

ygrid=4G. 

znin= "3. 

znax=3. 

ncontour=5. 


1 





Save the figure as a picture, leave the window settings as they are and superimpose the graph of z4. 




HlilH FiliD AUTD 3P 


This is a fine trick to show the obtain the whole solid of revolution. 
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(2) The other question might be more difficult: 



I'd like to evaluate the end points of intervals of convergence and I need the harmonic / alter- 
nating harmonic series. 

DNL: See first some sums of series calculated by the “Big Brother” DERIVE. Then I try to teach the 
Tl some rules: 

2 
7T 

n 

— co (-1) 

#4: 1 = - LNC2) 

n=l n 

co n co 

#5: T (-1) -n = Y n-COSOr-n) 

n=l n=l 


co 1 

I — 
#3 : n=l 3 

n 


2 

co n 

# 6 : I 

n-1 n + 2 



isum(u, v) 

Func 

Local u_, s_, d_, k_, d2 , res 
u / (1 /v)^u_: sign (u_)^s_ 
d(abs (u_) , v)^d2 
when (d2=0, 0,1 ,1)->d2 
string( ( "1 ) A v)-k_ 
stri ng (s_)^d_ 

If ( d_= " 1 " or d_= " - 1 " ) and d2=0 Then 
s_*oo^res : Goto end 
Endlf 

If d_=k_ and d2=0 Then 

u_*"I n (2) / ( "1 ) A n->res : Goto end 
Endlf 

when (I(u , v , 1 ,°°)=0 , 0 , 1(u , v , 1 ,°°) , "di v or ?")^res 

Lbl end 

res 

EndFunc 
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Heinrich Ludwig, Germany 

Lieber Herr Bohm, 

ich wiirde Sie gerne eine Kleinigkeit zu einem von DERIVES internen Algorithmen fragen. Vielleicht 
wissen Sie es und es kostet Sie nicht viel Miihe, mir einen Hinweis zu geben. 

Wie berechnet DERIVE die Wurzel einer ganzen Zahl? DERIVE muss einen sehr effizienten Algo- 
rithmus verwenden, um die quadratischen Faktoren einer ganzen Zahl zu erkennen. 

Multipliziert man z.B. a := 2 * MERSENNE(12) * MERSENNE(12) zu einer 77-stellige Zahl aus und 
bildet anschlieBend SQRT(a), so hat sogar mein Uralt-Rechner nach 0,02 s das Ergebnis gefunden. 
Eine Primfaktorzerlegung von a ist mit 0,1 1 s auch erstaunlich schnell, aber eben nicht schnell genug. 
DERIVE muss andere Wege kennen. 

Fur einen Hinweis, auch auf weiterfuhrende Literatur, ware ich Ihnen sehr dankbar. 

Mit freundlichen GriiBen 


Heinrich Ludwig 

Heinrich Ludwig asks which special algorithm is used by DERIVE for detecting quadratic factors of an 
integer. If you multiply a := 2 * MERSENNE(12) * MERSENNE(12) which results in a 77 digit number 
and then calculate SQRT(a) even Heinrich's very old PC needs only 0.02 sec to find the result, but if 
you factorize a then it needs 0.1 1 sec, which is also very fast, but significantly slower. DERIVE must 
know other methods. The DNL-reader knows that Johann Wisenbauer is the specialist for questions 
like these and he answered: 


Lieber Josef, 

Das ganze fallt unter das Thema 

detecting perfect powers. 

Herr Ludwig sollte also diese Stichworte in google eingeben und das Material dazu sichten, insbeson- 
derer auch die Arbeit 

• Daniel J. Bernstein (1998). " Detecting perfect powers in essentially linear time M . Mathematics 
of Computation 67 (223): 1253-1283. 

Wenn man jetzt von algorithmischen Feinheiten absieht, steckt nicht viel mehr dahinter, als dass man 
fur die in Frage kommenden n die nachstkleinere ganze Zahl zur n-ten Wurzel berechnet (meist mit 
Newton oder einer Modification davon) und dann schaut, ob die vorgelegte Zahl die n-te Potenz da- 
von ist. Wenn also im Zuge des Faktorisierungsprozesses einmal eine perfekte Potenz erhalten wird, 
wie im obigen Beispiel nach dem Heraussieben aller Primfaktoren bis 1024, dann geht der Rest sehr 
schnell. 

Herzliche GriiBe, 

Hannes 

This problem is part of the objective 

detecting perfect powers. 

Mr Ludwig should use Google and enter these key-words. He will find some materials, in particular the 
paper 


Daniel J. Bernstein (1998). " Detecting perfect powers in essentially linear time ". Mathematics 
of Computation 67 (223): 1253-1283. 
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Two-Stage Least Squares Linear Regression 

by MacDonald R. Phillips 


don.phillips@qmail.com 


The most widely used single-equation method for estimating simultaneous systems of equa- 
tions is two-stage least squares. Let Y be the endogenous or dependent variables in the 
system and X the exogenous or predetermined variables. The equations to be estimated are 
of the form: 

y=Y -(3 + X -y+u 
#5: 1 1 

where y is an n by 1 vector of observations on the "dependent" variable, Yjl is an n by g 
matrix of observations on the other endogenous variables included in the equation, (3 is the g 
by 1 vector of coefficients associated with Y|1 , X|1 is the n by k matrix of observations on 
the predetermined or instrumental variables appearing in the equation, y is the k by 1 vector 
of coefficients associated with Xjl , and u is the n by 1 vector of disturbances in the equa- 
tion. 


The problem of applying OLS to the above equation is that the variables in Yjl are corre- 
lated with u. The essence of two-stage least-squares regression is the replacement of Yjl 
by a computed matrix Y_hat|1 , where hopefully the stochastic element is purged, and then 
of performing an OLS regression of y on Y_hat|1 and Xjl . 

The matrix Y_hatj1 is computed in the first stage by regressing each variable in Yjl on all 
the instrumental variables in the complete model and replacing the actual observations on 
the Y variables by the corresponding regression values. Thus: 

-1 

#6: Y_hat = X-(X'-X) -X' -Y 

1 1 

where X = [Xjl X|2]. X is the n by K matrix of observations on aN the instrumental variables 
in the complete model, X|2 being the matrix of observations on those instrumental variables 
excluded from the equation under study. 


In the second stage y is regressed on 

Y_hatj1 and Xjl. 

The equation for the 2SLS esti- 

mates can be written as: 







-1 





-1 


Y '-X-(X'-X) -X'-Y 

Y l -X 




Y ’ -X-CX'-X) -X’-y 


1 1 

1 1 


' b ' 


1 

#7: 

X '-Y 

X '-X 

- 

. c . 

- 

X '-y 


1 1 

1 1 . 




1 


where b is the vector of coefficients on the other endogenous variables and c is the vector of 
coefficients on the predetermined variables in the equation, including the intercept. 
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The program TwoStage(eq_,endog_,exog_,data_,incept_:=1,vardef_:=1) computes the 2SLS 
regression coefficients, standard errors of the coefficients, t values and their probabilities, the 
ANOVA table, root MSE, R_square, and adjusted R_square. 

The inputs are: 

eq_: the equation to be solved for 

endog_: a vector of the endogenous variables 
exog_: a vector of the exogenous variables 

data_: the name of the dataset where the first row of the data matrix contains the names of 
the variables 

incept_: the default of 1 indicates that an intercept is to be computed for the equation; set 
incept:=0 if you do not want an intercept. 

vardef_: the default of 1 sets the variance denominator to the degrees of freedom; changing 
vardef_ to 0 sets the denominator to the number of observations. 


Note : According to SAS, if the no intercept option is set (incept_=0), the definition of the 
R_square statistic for two-stage least squares is changed to 1 - (Residal Sum of 
Squares/Uncorrected Total Sum of Squares). 

In addition, the program Model(mm_,endog_,exog_,data_,incept_:=1,vardef_:=1) sovles for 
all the of equations at once, mm is a matrix of the equations with one equation per row. 


Example 1: 

The following data are for a simplified model designed to explain variations in the consump- 
tion and prices of food. The data are from Kmenta, pp. 563-65. 

q = food consumption per head 
p = ratio of food prices to general consumer prices 
d = disposable income in constant prices 

f = ratio of preceding year's prices received by farmers to general consumer prices 
a = time in years 

The endogenous variables are q and p. The exogenous variables are d, f, and a. 

Estimate the following equations: 


(1) q = yO + pi * p + yl * d (the demand equation) 

(2) q = yO + pi * p + yl * f + y2 * a (the supply equation) 


Enter the data and then use the Model() program: 

q = p + d 
q = p + f + a 


#9: Model 


, [q, p], [d, f, a], data 
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#8: data := 


98.485 

100.323 

87.4 

98 

1 

99.187 

104.264 

97.6 

99.1 

2 

102.163 

103.435 

96.7 

99.1 

3 

101.504 

104.506 

98.2 

98.1 

4 

104.24 

98.001 

99.8 

110.8 

5 

103.243 

99.456 

100.5 

108.2 

6 

103.993 

101.066 

103.2 

105 .6 

7 

99.99 

104.763 

107.8 

109.8 

8 

100.35 

96.446 

96.6 

108.7 

9 

102.82 

91.228 

88.9 

100.6 

10 

95 .435 

93.085 

75 .1 

81 

11 

92.424 

98.801 

76.9 

68.6 

12 

94.535 

102.908 

84.6 

70.9 

13 

98.757 

98.756 

90.6 

81.4 

14 

105 .797 

95 .119 

103.1 

102.3 

15 

100.225 

98.451 

105.1 

105 

16 

103.522 

86.498 

96.4 

110.5 

17 

99.929 

104.016 

104.4 

92.5 

18 

105 .223 

105 .769 

110.7 

89.3 

19 

106.232 

113.49 

127.1 

93 

20 


MODEL: q 


'arameter 

Value 

Standard Error 

1(17) 

pi 

-0 , 2437076S9S 

0.09607652762 

2 , 536599686 

yo 

94. 614 85 3 36 

7 . 887 36 31 89 

11.99575207 

yi 

0 , 3143821013 

0.04674526445 

6 , 725432084 

Source 

DF SS 

MS 

F 


Prob(t) 


Model 2 202 . 7676225 101.3S3S112 26.4446926 

Residual 17 65. 17469563 3.S33S05625 

Corrected Total 19 267.94231S1 


Prob(F) 

DO 06 04 39446 5 8 


Standard Error R A 2 Adj_R A 2 

1.958010629 0.7567584841 0.728141835^ 


q = - 0. 2437076898- p + 0 . 3143821013 ■ d + 94.61485336 
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2 nd part of the output of simplification of ModelQ: 


# 10 : 


MODEL: q 


Parameter 

Value 

Standard Error 

t(16) 

Prob(t) 

PI 

0.240567722 

0.09950905774 

2.417545974 

0.02792589088 

*0 

49,44820072 

11.95947271 

4,134647227 

0.0007779536442 

*1 

0.2528436233 

0.09923148057 

2.548015248 

0.02148917651 

*2 

0.2560236611 

0.04704922378 

5.441612859 

0.00005432925863 


Source DF SS MS F Prob(F) 

Model 3 172 .1288502 57.37628342 9.581330838 0,0007369271092 

Residual 16 95.81346791 5.988341744 

Corrected Total IS 267.9423181 


Standard Error R A 2 Adj_R A 2 

2 .447108834 0,6424100957 0,5753615886 . 

q = 0. 240567722- p + 0. 2528436233- a + 0. 2560236611- f + 49.44820072 


Example 2: 

The next example is based on Klein's model 1 (1950). The endogenous variables are c, p, 
w, i, x, wsum, k, and y. The exogenous variables are klag, plag, xlag, wp, g, t, yr. 

yr = year - 1931 
c = consumption 
p = profits 

w = private wage bill 

i = investment 

x = private production 

wp = government wage bill 

g = government demand 

t = taxes 

k = capital stock 

wsum = total wage bill 

plag = profits lagged 

xlag = private product lagged 

klag = capital stock lagged 

y = c + i + g - 1 (national income) 

Estimate the model: c = f(p, plag, wsum), i = f(p, plag, klag), and w = f(x, xlag, yr). 
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Full matrix contains 21 rows: 


yr 

c 

P 

w 

i 

X 

wp 

9 

t 

k 

wsum 

plag 

xlag 

klag 

y 

-10 

41.9 

12.4 

25.5 

-0.2 

45.6 

2.7 

3.9 

7.7 

1B2.6 

28 . 2 

12.7 

44.9 

182.8 

37.9 

-9 

45 

16.9 

29.3 

1.9 

50.1 

2.9 

3.2 

3.9 

1B4.5 

32.2 

12.4 

45.6 

182.6 

46.2 

-8 

49.2 

18.4 

34.1 

5.2 

57.2 

2.9 

2 .8 

4.7 

1B9.7 

37 

16.9 

50.1 

184.5 

52.5 

-7 

50.6 

19 .4 

33.9 

3 

57.1 

3.1 

3.5 

3 .8 

152.7 

37 

IS. 4 

57.2 

189.7 

53.3 

-6 

52.6 

20.1 

35.4 

5.1 

61 

3.2 

3.3 

5.5 

157.8 

38.6 

15.4 

57.1 

152.7 

55.5 

-5 

55.1 

19 .6 

37.4 

5.6 

64 

3.3 

3.3 

7 

203.4 

40.7 

20.1 

61 

157.8 

57 

-4 

56.2 

15.8 

37.9 

4.2 

64.4 

3.6 

4 

6.7 

207.6 

41.5 

15.6 

64 

203.4 

57.7 

-3 

57.3 

21.1 

39.2 

3 

64.5 

3.7 

4.2 

4.2 

210.6 

42.9 

15.8 

64.4 

207.6 

60.3 

-2 

57. S 

21.7 

41.3 

5.1 

67 

4 

4.1 

4 

215.7 

45.3 

21.1 

64.5 

210.6 

63 

-1 

55 

15.6 

37.9 

1 

61.2 

4.2 

5.2 

7.7 

216.7 

42.1 

21.7 

67 

215.7 

53.5 

0 

50.9 

11.4 

34.5 

-3.4 

53.4 

4 .8 

5.9 

7.5 

213.3 

39.3 

15.6 

61.2 

216.7 

45.9 


Apply model() again: 



f 

c = p + plag + wsum 

'Y 

#12: Model 


i = p + plag + klag 

, [c, p, w, i, x, wsum, k, y], [klag, plag, xlag, wp, g, t, yr], klein 



w = x + xlag + yr 

J 


I'll show only the resulting model equations: 

[c = 0. 01730294013- p + 0 . 8101827436 .wsum + 0 . 2162334038- p lag + 16.55475199] 
[i = 0. 1502217351- p + 0 . 6159436535 - plag - 0 . 1577876492 - klag + 20.27821173] 
[w = 0.43885 8637 -x + 0 . 1466742262 -x lag + 0. 1303957913 -yr + 1.500299135] 
Same data with intercept option = 0: 




c = p+ plag + wsum 

> 

#14: Model 


i = p + plag + klag 

, [c, p, w, i , x, wsum, k, y], [klag, plag, xlag, wp, g, t, yr], klein, 0 



_ w=x+xlag+yr 



#19: [ c — 0.15 83375 706- p + 1 . 121162681 -wsum + 0 . 2651214661 - plag] 
#20: [i = 0.445965 7099- p + 0 . 3685120348- plag - 0 . 06148091416 - klag] 
#21: [w = 0.44445 845 25- x + 0 . 166299164 -x lag + 0 . 1129951702 -yr] 


References: 

J. Johnson, Econometric Methods, 2nd ed., McGraw-Hill Book Co., New York, 1972. 

Jan Kmenta, Elements of Econometrics, MacMillan Publishing Co., New York, 1971. 

L. Klein, Economic Fluctions in the United States: 1921-1940, John Wileys and Sons, New York, 
1950. 

J.K. Binkley and Nelson, G. (1984), "Impact of Alternative Degrees of Freedom Corrections in Two 
and Three Stage Least Squares," Journal of Econometrics, 24, 3, 223-233. 
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Heron for Square Roots 

Thomas Himmelbauer, Austria 

Thomas is one of the few teachers in Austria working with MuPAD in classroom. I am very grateful 
that he sent some MuPad papers with one lesson among them (Schulubung = school exercise): 


3 

Heronverfahrenl .mn* - MuPAD Pro 4.0 

Datei Bearbeiten 5uchen Ansicht Einfugen Formatieren Notebook Fenster Hilfe 

y 

3 sl Q & 

^ f 3 3 | ifc H> 0 @ 0^? 

Ge 


15. Schulubung 


H e ro nsch e F I ii ch e nfo rm e I : 


Diese Fortnel client zur nahrungsweisen Berechnung von Quadratwurzeln. 
Es handelt sich um ein rekursives Verfahren. 

This formula is used for approximative calculation of square roots, 
it is a recursive procedure. Demonstration of an example 

Erklarung anhand eines Beispiels: 
sqrt (51) 

' ^ 

1 . Schritt: 

Wir wahlen ein Rechteck, dessen Flache 51 groft ist. 

Step 1: Choose a rectangle with area = 51. 

~al:=l:bl:=51:al*bl 


51 


2. Schritt: 

Wir nehmen den Mittelwert der beiden Rechteckseiten als eine neue Rechteckseite, 
die andere bestimmen wir so dass der Flacheninhalt jeweils 51 ist. 

Step 2: We take the average of both sides as one new side and determine the second one 
in order to obtain again area = 51. 

~ a2 : = (al+bl) /2 :b2 : =51/a2 : a2 *b2 
51 


3.Schritt:| 

Wir nehmen den Mittelwert der beiden Rechteckseiten als eine neue Rechteckseite, 
die andere bestimmen wir so, dass der Flacheninhalt jeweils 51 ist. 

Step 3: Repeat step 2. 

~ a3 : = ( a2 +b 2 ) /2 :b3 : =51/a3 : a3 *b3 
51 

4. Schritt: 

Wir nehmen den Mittelwert der beiden Rechteckseiten als eine neue Rechteckseite, 
die andere bestimmen wir so, dass der Flacheninhalt jeweils 51 ist. 

Step 4: Repeat step 2. 

"a4:= ( a3 +b 3 ) /2 :b4 : =51/a4 : a4*b4 
51 

"float (b4 A 2) 

33.47S25S96 
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5. Schritt: 

Wiederhole Schritt 4. 

Step 4: Repeat step 2. 

~a.5:= (a4+to4) /2 : toS : =5 1/ a5 : a5 *1 d5 
51 

"float (b5 A 2 ) 

42 .80600931 

6. Schritt: 

Wiederhole Schritt 5. 

Step &: Repeat step 2. 

"a6: = (a5+to5) / 2 : to 6 : =5 1/ a6 : a6 6 
51 

"float ( a6 A 2 ) 

51.02465678 


Die Rechtecke nahern sich langsam einem Quadrat. 

Die Rechteckseiten nahern sich langsam der zu berechnenden Wurzel. 

The rectangles app>roach slowly a square and its sides approach the square root. 
plot (plot: : Rectangle (0 . . al, O . . to 1 ) , 


plot : 

: Rectangle (O . 

. a2 , 0 . 

■ b2) , 

plot : 

: Rectangle (O . 

. a3 , O . 

■D3) , 

plot : 

: Rectangle (O . 

. a4 , O . 

■b4) , 

plot : 

: Rectangle (0 . 

. a5 , 0 . 

■ bS) , 

plot : 

: Rectangle (□ . 

■ a6, □ . 

.b6) 


, Scaling=Constrained) 



X 


D 
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Da© Heronverfahren laftt sich durch folgendes Programm automatisieren! 

The procedure can be automated by the following program: 

Heron: =pr oc (A) 
begin 

M: ^matrix ( 7 , £ ) ; 

M[l,l] : —A; 

M[l,£] :=1; 

for i from £ to 7 do 

M[i, 1] : =float (l/£* (M [i-1, 1] +M [i-1, £ ] ) ) ; 

M [i, £ ] :=A/ float (l/£* (M [i-l, 1] +M [i-l, £ ] ) ) ; 
end_f or 
end_pr oc 

proc Heron(A) ... end 


Heron (£31) 

15. 19207849 


M 


/ 231 1 \ 

116.0 1.99137931 

58.99568966 3.915540294 
31.45561497 7.343680935 
19.39964795 11.90743257 
15.65354026 14.75704512 
\ 15.20529269 15.19207349 / 



X 


The following screen shots are from two MuPAD animations: 

The harmonic motion and “Sun - Moon and Earth in Motion” 
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Comments on demonstrating Heron’s Rule on the TI and with DERIVE. 

Performing the recursive procedure on the TI is very intuitive and informative as well in the Home 
Screen as applying the Sequence Mode. 


f » In 1 gebralca fclotherlprgn I die 1 ean Up! 1 

■51 ■+ a : 

1 ■* b 


1 

_ a + b _ 

■ 2 * a 

•4* b 


1.961538 

_ a + b _ 

' 2 * a 

= 4 ab 


3.647868 

- a + b + a 

= 4- 


5.78604G 

<a+b>/2+a : 51/a-Hbl 

MAIN 

RAD AUTD 5D 

4/50 



f » In 1 gebralca fclotherlprgn I die 1 ean Up! 1 

■ 51 * a : 1 ■* b 

1 

<a+b>/2+a:51 .ka+b 

MAIN RAD AUTD 5D 

i/50 



Not so easy, but of didactical value is defining the rule in the SEQUENCE Mode which unveils the 
recursive nature of the procedure in another way. We show how to approximate V231 . 


PJSJsetupT::?! ( 

- FH Y 

Header h; 

Al. 

™V>-T 1 

n 

Ul 

u2 




0.00000 

231.000 

1.00000 




1.00000 

116.000 

1.99138 




2.00000 

58.9957 

3.91554 




3.00000 

31.4556 

7.34368 




4.00000 

19.3996 

1 1 . 9074 




5.00000 

15.6535 

14.7570 




6.00000 

15.2053 

15. 1921 




7.00000 

15. 1987 




u2 <n >= 

=15.198682717465 



MAIN 

RAD AUTD 

SEG 



In Derive we can not only realise the recursive nature by applying the ITERATES-command, but also 
animate the approximating process using Peter Schofield's Slide Show tool. (For details about the 
Slide Show see the following contribution.) 


# 1 : 


heronfz , 


n) 


v + v 

1 2 


:= ITERATES 


2 


v 


2 ■ z 


V + V 

1 2 _ 




v, [z, 1], n 


J 


heron(Sl, 6) 


SI 

1 

26 

1.961538461 

13. 9S0 76923 

3.647S6795 

S. S1431SS9 

S . 786040007 

7 . 30017929S 

6 . 9S61297S 

7 . 1431S4S39 

7.139702734 

_ 7 . 141428637 

ysi 

7 . 14142S219 


# 5 : 
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#6 creates the family of rectangles which converge to a square. 


#6: heron_sq(z, n) := ITERATES 


r 


0 


0 





V 

2,1 

+ V 

3. 

2 

0 






2 



" 0 

0 " 



V 

2,1 

+ V 

3. 

2 

2 ■ z 

, V, 

z 

0 

1 

, n 



2 


V + V 

2,1 3,2 


£- 

_ 0 

_L 

1 _ 




0 


2 ■ z 







V + V 

2,1 3,2 . 



j 


Slide_Show(v, k) := VECTOR (v -IFCk = i, 1, ?) , i, 1, DIMM) 
#7: i 

Slides_ShowM k) := VECTOR (v -IFCk > i, 1, ?) , i, 1, DIMM) 

#8: i 

#9: root51 := heron_sq(51, 10) 

# 10 : 


#6 creates the family of rectangles which converge to a square. Highlight #10, switch to the 2D-Plot 
Window, create a slider bar for k with 1 < k < 1 1 with 10 intervalls. (Points small and connected) 


1.00 


-> 


11.00 


12345678 9 10111213141516171819 20 2122 23 24 25 26 27 



1 2 3 4 5 6 7 8 9 10111213141516171819 20 2122 23 24 


Now you can stepwise follow how the initial rec- 
tangle changes into a square. (You could add the 

shaded square with sides Vil in order to make the 
approximating procedure more impressive. 

(On the right picture we are on the way to y/l3l .) 



1 2 3 4 5 6 7 8 9 1011121314151617181920 2122 2324 25 26 
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A Slide Show for the Slider Bar 

Peter Schofield, UK 


We had a seminar in Amstetten, Lower Austria, for evaluating several learning paths 
and discussing possible ways of integrating various technologies and media in math 
teaching and learning. During a break - breaks are a very important part of seminars 
for exchanging ideas and for receiving inspirations - one colleague asked for advice 
how to animate polygons inscribed in a circle to visualise approximation of n. She 
set up a VECTOR and could easily produce the polygons and plot one after the 
other. 


poly(n) := VECTOR 


f 


^ 2 «7T« 1 ^ 


r 2 «7T« 1 ^ ‘ 

> 


cos 


, SIN 


i I, 0 , n 






^ n ^ 


k. n j. 



polys(m) := VECTOR (poly (n) , n, 3, m) 


She wanted to install a slider bar for showing one polygon after the other. Each polygon appears as 
an element of the list polys(m). We both were not able to produce the performance of the polygons. I 
remembered a similar problem to show points of a list one after the other (DNL#62) but I couldn't 
make David Sjostrand's trick working (and didn't try Peter Schofield's idea ...). I knew that there are 
often problems occuring when applying the slider bar for a discrete variable. I wrote to Peter and de- 
cribed the problem. Peter came back to me after a couple of days - and here is his terrific slide show 
tool for two kinds of demonstrations. Josef 


SlideShow.dfw: a Derive 6 file for 2D-plotting a slide-show of regular n-sided polygons lead- 
ing to an approximation for a circle. March 2007. 

Peter Schofield, Trinity & All Saints, LEEDS - email: p.schofield@leedstrinity.ac.uk 


Si ide_Show(v , k) := VECTOR (v .IF(k = i, ?) , i, 1, DIM(v)) 

#1: i 

Slides_Show(v, k) := VECTOR (v -IF(k > i, ?), i, 1, DIM(v)) 
#2: i 


#3: poly(n) := VECTOR 


r 


r 2 -7T' 1 ^ 


r 2 ■ rr ■ i ^ 

> 


COS 


, SIN 


■ i i 0, n 

j ■ 





^ n J 


^ n j . 

j 


#4: polys(m) := VECTOR (poly(n) , n, 3, m) 

0 


#5: polytri(n) := VECTOR 


COS 


r 2 -rr - 1 ^ 

SIN 

^ 2 ■77* i ^ 



^ n j 


^ n j 


COS 


2‘TT- (i + 1) 


SIN 


2 -jr- (i + 1) 


i , 0, n - 1 


#6: polytris(m) := VECTOR(polytri(n) , n, 3, m) 


siide_show shows one polygon after the other and siides_show superimposes one graph after 
the other. 

Select 2D-plot Window (Approx Before Plotting (ON)). 

Insert>Sl ider-Bar> 3<k<36 with 33 intervals. 

2D-plot and slide the following: 
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#7: Si nde_Show(polys(36) , k - 2) 

#8: Si ndes_Show(polytns(36) , k - 2) 

#9: Si "ides_Show(polys06) . k - 2) 

If you don't believe me then you are friendly invited to make a try - you can download all 
files - this works wonderful. (Don't forget to set Approximate before Plotting on and to have 
the points connected, point size small is recommended.) 


Peter closed his email by pointing out that this tool might find other applications, too. I took 
his words for serious and accomplished the inscribed polygon file with circumscribed poly- 
gons and a visualisation of the approximation process for n. 

In- and circumscribed polygons with variable radius r 

r- cos 


#10: poly_-in(r, n) := VECTOR 



# 11 : 


poly_out(r , 





#12: x + y =9 

#13: poly_ins := VECTOR (poly_in (3 , n) , n, 3, 36) 

#14: Si ide_ShowCpoly_ins , k - 2) 

#15 : poly_outs := VECTOR (poly_out (3 , n) , n, 3, 36) 

#16: Si ide_Show(poly_outs , k - 2) 

#1 7 and #1 8 plot segments of height approximating tt - converging from below and from above 
It is recommended to plot #14 and #1 7 in the same colour and #1 6, #1 8 in another one. 

#1 9 gives a horizontal line y = u 

' 5 o 

#17: pi_in(k) := k f 2-jt ''j 



#19: 7T 
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Here are some screen shots: 




; J f 7 


k = 11.00 \Xj\ 



n 



When I read Peter's note I immediately had the idea to use the slide show for demonstrating 
integration. More and more teaching software tools have animated Riemann sums imple- 
mented - we can now implement the animation by ourselves - thanks Peter's slide show 
tool. 


Another application: 

Left and right Riemann sums 


#20: leftsums(u, v, a, b, n) := VECTOR 


VECTOR 


(i - 1Mb - a) 


(i - 1Mb - a) 
n_ 

Mb - a) 


SUBST u, v, a + 


SUBST u, v, a + 


n_ 


(i - 1Mb - a) 
n_ 

(i - 1) • (b - a) 
n_ 


Mb - a) 


#21: Is := leftsums 


+ 1, x, 1, 5, 50 


# 22 : 


|S~Hde_Show(1s, k)| 
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Insert a slider bar for k = 1 ... 50 with 49 intervals, plot the function and then plot #22. 



Do the same for thr right sum (or the midpoint sum or the trapezium sum) and visualise the approach 
from both sides to the area. 



1 2 3 4 5 6 7 S 9 10 11 12 
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Black and White 
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Impressions in Pastel 
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3D Heat Waves 






F(xy> - |Kt‘(«/3> (COS<y) - SIM<y)) 
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Titbits from Algebra and Number Theory (33) 

by Johann Wiesenbauer, Vienna 


This column deals with permutations on a finite set, say {l,...,n}, which have received a 
lot of attention, as they occur in many places in mathematics and elsewhere. Just as a 
starter let me quote the following nice problem by Ivors Peterson at 
http://my.opera.com/mathmhb/bloq/show.dml/415087 

"The names of 100 prisoners are placed in 100 wooden boxes, one name to a box, and the 
boxes are lined up on a table in a room. One by one, the prisoners are led into the room. 
They may look into up to 50 of the boxes to try to find their own name, but must leave 
the room exactly as it was.The prisoners are permitted no further communication after 
leaving the room. They do have a chance to plot a strategy in advance. Good thing. 
Unless they all find their own names, they will all be executed. If each prisoner exam- 
ines 50 boxes at random, the probability of the group's survival is a miniscule 1/2^100 
[...]. However, there is a strategy that the prisoners can use to i ncrease the probability 
of success to more than 30 percent. Incredible but true! The trick is to find this 
strategy." 

Believe me, the solution is a real jawdropper! Hence, before reading on, it might be a 
good idea to stop here and have a try at this nice problem yourself! 

As you may imagine by now, it has to do with special properties of permutations. First 
the prisoners have to agree on some bisection between the boxes and their names. For 
the sake of convenience, we assume that both the boxes and the prisoners are num- 
bered with the numbers 1,2,. .,100. Then each prisoner, if it is his turn, goes to the box 
whose number corresponds to his number, reads the name contained in it and interprets 
it as a pointer to the next box, again reads the name contained in it and uses it as a 
pointer to the next box etc. 

It should be clear by now that this strategy is successful if and only if the underlying 
permutation (consisting exactly of the pointers above!) has got no cycles of length 
greater than 50. For people, who are good at combinatorics, in particular as regards the 
manipulation of generating functions, it is not too hard to identify the probability that a 
permutation on {l,2,...,n} has no cycle longer than k as the coefficient of x n in the taylor 
series of exp(x + x 2 /2+...+x k /k). 

This is where Derive comes into play! Again, I invite you to write a routine for the com- 
putation of this probability for general values of n and k yourself, before looking at my 
solution below. (It is the final solution after quite a few attempts before had ended in 
an overflow of memory or took too much time!) 


probO(n, k, n_ := 0, r_ := 1, s_ := 0, t_) 
Prog 

t_ := rCx A k_/k_, k_, 1, k) 

Loop 

s :+ P0LY_C0EFF(r_, x, n) 

n_ : + 1 
If n_ > n 


RETURN s 

r := REMAINDER(r_-t_, x A (n + l))/n_ 
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TABLE(probO(4, k) , k, 1, 4)’ 


Can you check this 
manually? 


probOCLOO, 50) = 


21740752665556690246055199895649405434183 


69720375229712477164533808935312303556800 


(3.75s) 


( 2174075 2665 5 5669024605 5199S9 5 6494 054341 S3 ") 

APPROX 100 = 31.18278206 

^ 69720375229712477164533808935312303556800 j 

If all this looked a little frightening to you, here comes a bit of relief. Provided that 
2k + l>n , which is certainly fulfilled in the case at issue, there is a much simpler for- 
mula to compute those probabilities. For the derivation of this formula, lets consider 
first the probability that a given permutation on {1,2, „,n) has a cycle of length j with 
j > n/2 . Since there are (j — 1) ! cycles of length j (note that we may always assume 
w.I.o.g. that the smallest element of the cycle is also the very first one) and 
n(n-l)...(n-(n-j) + l) is the number of possible assignments for the remaining elements, 


we get 


n (n-l)...(j + l)(j-l)! 


for our probability, which simplifies after cancelling to the surprisingly simple expres- 
sion 1/j. Hence, by summing up these probabilities for j = k + l,...,n and forming the 
counter probability, we get the value 

1- V -»l-ln- 

j±r. j k 

for the probability that there are no cycles of length greater than k for any given k 
with 2k + l>n. (Keep always in mind the latter restriction for the validy of this for- 
mula!) If you know the approximation In 2 « 0.693 (and every educated mathematician is 
expected to know this value by heart, isn't he?), then you get the rather good - though 
a little too small - approximation 0.307 for the probability we are aiming at, where 
n/k = 2 . 

Using this knowledge, we can implement a far better version of probO(n,k). (Well, as its 
name already indicates, we had something better in mind all the time!) 

prob(n, k) := 

If 2-k + 1 < n 
prob0(n, k) 

1 - j_, k + 1, n) 


prob(100, 50) 


2174O75266555669O246O55199S956494054341S3 


69720375 2297124771645 33SOS9353123O3556SO0 


(0.000s !!) 


This shows again that what we really need is not faster computers, but better algo- 
rithms! 
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Ok, let's turn to a different topic, which is the efficient traversing of all permutations 
in a given order (in the following always the lexicographical order). Since permutations 
on {l,2,..,n} have the annoying property that their total number n! increases very much 
with n, it is impossible for bigger values of n, say n > 9, to have a copy of each one in 
memory simultaneously. For most applications it suffices though, to find the "next" 
permutation of a given permutation wrt to the given order. Furthermore, usually certain 
conditions must be fulfilled, i.e. not all numbers are allowed in all positions. (Note that in 
the following routines the list of "forbidden" elements for each position can be found in 
the vector v.) Hence, we have to deal with the problem to find the next "valid" permuta- 
tion. 

It turns out though that we first need another important routine that will find for us 
the very first valid permutation under the given conditions and wrt the given order. 

f i rst_perm(s , v := [], f_, n_, s_, v_) := 

Prog 

If s = [] 

RETURN [] 

If v = [] 

v := VECTOR ([], k_, DIM(s)) 

f_ := SELECT (-. MEMBER?(f , FIRST(v)), f_, s) 

If f- = [] 

RETURN ? 
f_ FIRST (f_) 
v_ := REST (v) 

s_ := SELECTfs * f_, s_, s) 

If f i rst_perm(s_, v_) = ? 

first_perm(s, AD30IN(AD30IN(f_, FIRST(v)), REST(v))) 

AD30IN(f_, f "i rst_perm(s_, v_)) 

Note that if v is empty (a rather sloppy way to say that v consists of empty lists only, 
i.e. no elements are forbidden in any position), then the very first permutation is always 
the identity mapping on (1,2,... ,n). 

As for next_perm(p,v), assuming that p is a valid permutations wrt v, we are looking for 
the rightmost position of a possible transposition. Let i and j with i<j be the correspond- 
ing positions in p. Then we first perform that transposition and then look for the very 
permutation that leaves the positions 1,2,..,I unchanged. For the rest, I leave it to you 
to find out the details. 


v := 

Prog 

If v = [] 

1 

cr 

1 1 

1 1 

1 

Cl 

v := VECTOF 

?([], k_, DIH(p): 

p := REVERSECf 


v := REVERSED 
Loop 

0 


p_ := AD30IN(FIRST(p), p_) 
p := REST (p) 

If P = [] 

RETURN ? 

v_ := AD30IN(FIRST(v), v_) 
v := REST (v) 

q_ := SELECT (FI RST(p) < q_ a - MEMBER?(q , FIRST(v)), q_, SORT(p_)) 

Loop 

r_ AD30IN(FIRST(p) l SELECT(p_ * FIRST(q_) , p_, p_)) 
r_ := f i rst_perm(SORT(r_) , v_) 

If r_ t ? exit 
q_ := REST (q_) 

If q_ = [] exit 
if q_ a [] 

RETURN APPEND(REVERSE(AD30IN(FIRST(q_) , REST(p))), r_) 
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ITERATE (next_perm(p_), p_, [1 7], 7! - 1) = [7, 6, 5, 4, 5, 2, 1] 

Note also the totally "index-free" style of programming in the programs above, some- 
thing I'm really very proud of! In my experience, people coming from other programming 
languages like Pascal, VB, and diverse C-dialects etc. have often really a problem to get 
used to this style. Ok, here we go again, you may be thinking! You are right, I presume 
that you are eager by now to see these routines at work. 

Hence, let's use them for example to classify all quasigroups (or latin squares, if you 
prefer that) on the set {l,2,...,n}. They are essentially (nxn)- matrices with the property 
that all rows and columns are permutations of {l,2,..,n}. In other words, each new row 
has the property that it is a permutation where the vector v of "forbidden" elements is 
exactly the transpose of the matrix built so far. The boolean variable r, which is set to 
true by default, determines whether we are only interested in so-called "reduced" qua- 
sigroups, which have the property that the first row and the first column are both the 
identity mapping. (Note that you can get all other quasigroups by first permuting the 
column and afterwards the rows, leaving the very first row fixed. In particular, the to- 
tal number of quasigroups is bigger by a factor n!(n-l)!.) 

quasigroups(n , r := true, s := [[]], p_, q_) := 

Prog 

If r a s = [[]] 

s := [[[1 n]]] 

Loop 

q_ := FIRST (s) 

If DIM(q_) = n 
RETURN s 
s := REST (s) 

p_ := f-j rst_perm([l n], q_’) 

Loop 

If p_ = ? exit 

If r a FIRST (p_) > DIM(q_) + 1 exit 
s := APPEND(s , [APPEND(q_, [p_])]) 
p_ := next_perm(p_, q_’) 

T 1 2 3 4 1 T 1 2 3 4 1 I" 1 2 3 4 1 T 1 2 34 

2 14 3 2143 2341 2413 

quasi groups (4) = , , , 

3412 3421 3412 3142 

_L 4 3 2 1 J 1.4 31 2 j [4 1 2 3 J 1.4 3 21 

DIM (quasi groups (4 , false)) = 576 
DIM(quasigroups(5)) = 56 

As the notation "quasigroups" already indicates, there is a close relationship of qua- 
sigroups to groups. More precisely, every group G is a quasigroup, where additionally the 
so-called associativity law (xy)z = x(yz) for all x,y,z e G holds. This then implies the 

existence of an identity element eeG such that ex = xe = x for all xeG and the exis- 

tence of an inverse element x'eG for each xeG such that x'x = xx' = e . (Try to prove 
this!) 
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Note that all 4-element reduced quasigroups above are actually groups. This is no longer 
true for 5-element reduced quasigroups (or slightly more generally, quasigroups with an 
identity element, which are also called loops), as you may check yourself referring to 
the complete list of those quasigroups below, which was established using a special aux- 
iliary routine for the output of such lists in matrix form. 

matrixCu, k) := 

Prog 

Loop 

If MOD(DIMCu), k) = 0 exit 
u := APPENDfu, [ m1 ]) 

VECT0R(u| [k ■ j_ + 1 k ■ j_ + k], j_, 0, FL00R(DIM(u) , k) - 1) 

matrix (quasi groups (5) , 8) 



But there are other well-known mathematical objects, which can also be viewed as spe- 
cial quasigroups, and this time they are even known to the proverbial man in the street. 
Just in case, that didn't ring a bell yet, I'm talking about sudokus. 

Well, as you may remember, we have already dealt with sudokus in this series, but using 
a completely different approach. Now we are viewing them as special quasigroups on the 
set {1,2,. „n}, where n must be of the form n = m 2 for positive integer m, and try to 
take advantage of our routines dealing with permutations. 
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As a rule of thumb, the run time is several minutes rather than several seconds (the 
example above took about 6 min on my PC). Hence, this is not meant as a sudoku-solver, 
but rather as a demonstration that our ideas work. On the other hand, it will solve any 
suduko reliably and also tell you all solutions, if there is more than one. (At least, this is 
the default setting, which be changed though by setting the parameter all to false.) 
Well, and not to forget, it also works for other values of m than the “classical" case 
m=3, at least in principle. Let's conclude with the following example, where m=4. 


sudokus(4, [], false) 


1 

2 

3 

4 

5 

6 

7 

8 

5 

6 

7 

8 

1 

2 

3 

4 

9 

10 

11 

12 

13 

14 

15 

16 

13 

14 

15 

16 

9 

10 

11 

12 

2 

1 

4 

3 

6 

5 

8 

7 

6 

5 

8 

7 

2 

1 

4 

3 

10 

9 

12 

11 

14 

13 

16 

15 

14 

13 

16 

15 

10 

9 

12 

11 

3 

4 

1 

2 

7 

8 

5 

6 

7 

8 

5 

6 

3 

4 

1 

2 

11 

12 

9 

10 

15 

16 

13 

14 

15 

16 

13 

14 

11 

12 

9 

10 

4 

3 

2 

1 

8 

7 

6 

5 

8 

7 

6 

5 

4 

3 

2 

1 

12 

11 

10 

9 

16 

15 

14 

13 

16 

15 

14 

13 

12 

11 

10 

9 


9 

10 

11 

12 

13 

14 

15 

16 

13 

14 

15 

16 

9 

10 

11 

12 

1 

2 

3 

4 

5 

6 

7 

8 

5 

6 

7 

8 

1 

2 

3 

4 

10 

9 

12 

11 

14 

13 

16 

15 

14 

13 

16 

15 

10 

9 

12 

11 

2 

1 

4 

3 

6 

5 

8 

7 

6 

5 

8 

7 

2 

1 

4 

3 

11 

12 

9 

10 

15 

16 

13 

14 

15 

16 

13 

14 

11 

12 

9 

10 

3 

4 

1 

2 

7 

8 

5 

6 

7 

8 

5 

6 

3 

4 

1 

2 

12 

11 

10 

9 

16 

15 

14 

13 

16 

15 

14 

13 

12 

11 

10 

9 

4 

3 

2 

1 

8 

7 

6 

5 

8 

7 

6 

5 

4 

3 

2 

1 
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A Conjecture - A Theorem? 


Heinrich Ludwig, Germam 


Lieber Herr Bohm, 

ich danke Ihnen fur Ihre schnelle Auskunft! Die Arbeit von Daniel J. Bernstein konnte ich liber das Internet 
beziehen. Mit dem Studium derselben muss ich mich leider noch gedulden. Meine Arbeit im Dienst fordert mich 
zur Zeit deutlich. 


Vielleicht interessiert es Sie, dass ich mit DERIVES Hilfe auf eine (mich) verbluffende Vermutung gestoBen bin. 
Um den Aufgabenfundus fur meinen Leistungskurs zu erweitem, habe ich Kreise im 2-dimensionalen Raum 
(wie iiblich in der Schule mit euklidischer Metrik) gesucht, die durch 3 Punkte mit Ganzzahl-Koordinaten gehen 
und die zugleich einen rationalen Radius haben. Man findet eine ganze Menge Losungen. Sortiert man sie der 
GroBe nach, so ergibt sich eine Folge von Briichen, die einigermaBen irregular erscheint. Neugierig geworden, 
habe ich DERIVE dann entsprechende Losungen im 3-, 4-, 5- und 6-dimensionalen euklidischen Raum suchen 
lassen. Das Erstaunliche ist, dass im IR_5 (aber noch nicht im IR_4!) eine gewisse Vollstandigkeit erreicht wird, 
die folgender Satz ausdriickt: 

Sei q eine naturliche Zahl und p’ := (q A 2)/2 falls q ungerade, bzw. p' := (q A 2)/4 falls q gerade, dann gibt es fur 
jede naturliche Zahl p > p' drei Punkte im IR_5, so dass der Kreis durch diese Punkte den Radius p/q hat. 

Zu einem Beweis fehlt mir leider jeder Anhaltspunkt. Trotz tagelanger Rechenzeit hat sich aber auch kein Ge- 
genbeispiel eingestellt. - Die Umkehrung des Satzes gilt nicht, es gibt auBer den eben beschriebenen Umkreisra- 
dien einige sporadische Werte, die sich einer naheliegenden Einordnung entziehen. Neugierig geworden, was 
hohere Dimensionen noch zu bieten hatten, untersuchte ich den IR_6 ebenso. Erstaunlich ist, dass im IR_6 keine 
weiteren Losungen hinzukommen. Hoherdimensionale Raume als den IR_6 sind schwer ausreichend zu untersu- 
chen, weil die Rechenzeit mit der Dimension dramatisch wachst. 

Ich raume ehrlichkeitshalber ein, dass ich nach meinen ersten Berechnungen auf Programme uberging, die in 
Delphi (Pascal- Abkommling) geschrieben sind, weil damit eine gegeniiber DERIVE 800-fach hohere Rechenge- 
schwindigkeit moglich war. Aber: DERIVE hat die entscheidenden Hinweise gegeben, die mich motiviert haben 
tiefer nachzubohren! 


Mit freundlichem GruB 


Ihr Heinrich Ludwig 


P.S. Meine ersten Ergebnisse habe ich (diese alle noch mit DERIVE ermittelt) in der bekannten ’’On-Line En- 
cyclopedia of Integer Sequences” von Neil Sloane unter den Indizes A 128006 ff. veroffentlicht. 


I try to translate Heinrich's interesting findings (Josef). 

I found a (for me) surprising conjecture (supported by DERIVE). Preparing problems for my classes I 
looked for circles having a rational radius and passing points with integer coordinates. One can find a 
lot of solutions. Sorting them in increasing order one obtains a sequence or fractions which seems to 
be quite irregular. Extending the research for 3, 4, 5 and 6 dimensional space I noticed that in IR_5 a 
certain completeness seems to be reached which can be expressed as follows: 

Let q a natural number and p’ :=q 2 /2 if q odd and p’=q 2 /4 if q even, then for every natural number p > 
p’ exist three points with integer coordinates in IR_5 such that the circle passing these points has a 
rational radius p/q. 

Unfortunately I don't have any idea for a proof. I worked for days and couldn't find any counter exam- 
ple. Reversion of the theorem is not valid. Then I extended my search to IR_6. It is surprising that 
there are no additional solutions appearing. Investigations of higher dimensioned spaces are difficult, 
because calculation times will increase dramatically. 

To be honest, I must admit, that after my first calculations I changed to programs written in Delphi, 
because calculation time is 800 times faster. But it was DERIVE which gave important hints, which 
inspired me to get deeper into the matter. 

I published my first results (all of them found with DERIVE) in the well known “On-Line-Encyclopedia 
of Integer Sequences" of Neil Sloane under the index numbers A1 28006 ff. 






